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- - ■ Abstract 

t^ ■ 

^\ I In this chapter, we consider the radiation reaction to the motion of a point-hke 

particle of mass m and specific spin S travehng on a curved background. Assuming 
, S = 0{Gm) and Gm <^ L where L is the length scale of the background curvature, 

(]j ' we divide the spacetime into two regions; the external region where the metric is ap- 

Q ■ proximated by the background metric plus perturbations due to a point-like particle 



and the internal region where the metric is approximated by that of a black hole 
plus perturbations due to the tidal effect of the background curvature, and use the 
technique of the matched asymptotic expansion to construct an approximate metric 
which is valid over the entire region. In this way, we avoid the divergent self-gravity 



(N 



> 

\^ ' at the position of the particle and derive the equations of motion from the consis- 

J;^ ■ tency condition of the matching. The matching is done to the order necessary to 

^^ I include the effect of radiation reaction of 0{Gm) with respect to the background 

metric as well as the effect of spin-induced force. The reaction term of 0{Gm) is 

t^^ ■ found to be completely due to tails of radiation, that is, due to curvature scattering 

^^ \ of gravitational waves. In other words, the reaction force is found to depend on the 
O 



entire history of the particle trajectory. Defining a regularized metric which consists 

of the background metric plus the tail part of the perturbed metric, we find the equa- 

t^, tions of motion reduce to the geodesic equation on this regularized metric, except 

for the spin-induced force which is locally expressed in terms of the curvature and 

. spin tensors. Some implications of the result and future issues are briefly discussed. 

rS 

C^ ■ §1. Introduction 

The problem of radiation reaction has long been one of the fundamental theo- 
retical issues in general relativity. Starting from the historical works of Eddington 
in his 1922 bookEP, Chandrasekhar and Espositon^ discussed the radiation reaction 
of the self-gravitating fluid emphasizing the importance of the time asymmetric part 
of the metric appearing in the post-Newtonian expansion, and Burke and Thorne™ 
found that the leading contribution from the time asymmetric part can be compactly 
expressed in the form of a resistive potential. 

The previous studies of radiation reaction B'u) were done under the assumption 
that the post-Newtonian expansion is valid. Here we consider this problem in the 
framework of linear perturbation theory in a general spacetime. A part of moti- 
vation is to give a rigid foundation of the method to solve the Einstein equations 
perturbatively as an expansion with respect to the perturbation caused by a point- 
like particle. Usually one adopts a point-like particle to represent a black hole or 



neutron star in the linear perturbation studies as was done in Chapter 1. Then one 
may pose several questions: Since the perturbed field diverges at the location of the 
point-like particle, is the approximation scheme of linear perturbation still valid? 
Does the point-like particle really represent a black hole or a neutron star? If it 
represents a black hole, the center of it is inside the event horizon, and then in what 
sense does 'the motion of the particle' make sense? Here we are going to clarify the 
meaning of the particle trajectory and derive the equations of motion including the 
effect of radiation reaction to the first non-trivial order. 

Before starting the discussion of the gravitational radiation reaction, it is worth- 
while to refer to the electromagnetic case in a fixed curved background spacetime 
which was discussed by DeWitt and Brehmea^ In the electromagnetic case, the total 
energy momentum tensor composed of the particle and field contributions satisfies 
the conservation law. The conservation law is integrated over the interior of a world 
tube with an infinitesimal length surrounding the particle orbit. The part of the 
integration which does not vanish in the limit of small tube radius is transformed 
into the surface integrations over both ends of the tube and over the surface of the 
tube by using the Gauss theorem. The integrations over the top and bottom of 
the tube, respectively, give the definition of the particle momenta at both ends and 
the difference between them represents the change of the momentum during this 
infinitesimal time interval, which is to be equated with the momentum fiow given by 
the integration over the surface of the tube. In this way the equations of motion are 
obtained. 

In the case of gravitational radiation reaction, it is possible to construct a con- 
served rank-two tensor defined ori the background spacetime, composed of the matter 
field and the metric perturbation^'. However, there is an essential difference between 
the electromagnetic and gravitational cases. In electromagnetism, we can consider 
an extended charge distribution which is supported by a certain force other than 
the electromagnetic field. Thus it is possible to assume that the charge and mass 
distributions of a point-like particle are not distorted by the effect of the radiation re- 
action. Therefore one may consistently assume that the momentum and the electric 
current of the particle are proportional to the 4-velocity of the particle. Moreover 
the electromagnetic charge e is not directly related to the energy momentum of the 
particle which is proportional to the mass m. Hence, even if the limit of zero particle 
radius is taken, the divergent self-energy (oc e"^) can be renormalized into the mass. 
In the case of gravitational radiation reaction, it is not possible to consider such an 
ideal point-like particle because every force field universally couples with gravity. 
Even worse, the role of e in electromagnetism is also attributed to m. Thus a simple 
renormalization scheme does not make any sense. 

In order to deal with the gravitational case, we use the matched asymptotic ex- 
pansion technique that has been studied by many authors (e.g., D'EathQ) and Thorne 
and HartleQ') in the context of the post-Minkowski (or post-Newtonian) approxima- 
tion. We assume that the metric sufficiently far from the particle is approximated 
by the perturbation on the background spacetime generated by a point-like particle. 
We call this the external metric. We also assume that the internal metric which 
describes the geometry around the particle is represented by a black hole metric of 
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small mass in the lowest order approximation. As the particle moves in the curved 
background, the internal metric suffers from the tidal distortion. Thus both internal 
and external metrics are constructed perturbatively. The expansion parameters for 
the internal and external metrics are, however, different. We call this construction of 
the metric in the internal region the internal scheme and that in the external region 
the external scheme. Assuming the existence of the matching region where both 
schemes are valid, the internal and external metrics are expanded there as double 
series with respect to the two expansion parameters. Then the terms in these series 
are labeled by two indices which denote the powers of the two expansion parameters. 
Equating them order by order, we obtain the matching condition, through which one 
scheme determines the boundary condition of the other and vice versa. 

Using the matched asymptotic expansion to the first non-trivial orders of the 
expansion parameters, we present two different derivations of the equations of motion 
with the radiation reaction force of 0{Gm); (1) by means of an explicit construction 
of the metric, and (2) by using the so-called laws of motion and precessionLP. 

As mentioned above, in constructing the internal metric, the tidal distortion of 
the geometry is taken into account by the perturbation of the black hole. In the 
method (1), we set the gauge condition in the internal metric so that J = and 
1 linear homogeneous perturbations of the black hole vanish since they are purely 
gauge degrees of freedom as long as both the mass and angular momentum of the 
black hole stay constant. Applying a limited class of coordinate transformations 
that keep the meaning of the center of the particle unambiguous, the internal metric 
is matched with the external one in the matching region. Then we find that for 
a given trajectory of the particle a consistent coordinate transformation does not 
always exist, and this consistency condition determines the equations of motion. 

In the method (2), not all the metric components are evaluated in both schemes 
independently but we assume the existence of a coordinate transformation that gives 
a relation between the internal and external metrics. Once we know some metric 
components in one scheme, the counter parts in the other scheme are obtained from 
the matching condition. At this stage, the gauge condition is not fixed in a unique 
manner. The coordinate transformation between the internal metric and the external 
metric is chosen so that some of the metric components that are evaluated in both 
schemes are correctly matched in the matching region. Substituting the metric con- 
structed in this way into the Einstein equations, we obtain the consistency condition. 
There is a convenient method to extract out the information about the equations 
of motion from the consistency condition. Namely to use the laws of motion and 
precession introduced by Thorne and HartleQ^ . The laws of motion and precession 
are derived from the non-covariant but conserved form of the Einstein equations. 

The resulting equations obtained from both derivations are the same, although 
the strategies are quite different. In the method (1), the metrics in both schemes are 
calculated independently by using the Einstein equations. The matching condition is 
used to obtain the consistency conditions, which in turn give the equations of motion. 
On the other hand, in the method (2), the matching condition is used to construct 
the metric. The consistency condition is derived by requiring that thus obtained 
metric satisfies the Einstein equations. The meaning of the matching condition in 



deriving the equations of motion is clearer in the method (1) than in (2), but the 
method (2) is much simpler and straightforward than the method (1) as we shall see 
in the following. 

The organization of this chapter is as follows. We use the terminology 'a 
monopole (spinning) particle' to refer to a particle which represents a Schwarzschild 
(Kerr) black hole. In section y, the matched asymptotic expansion technique is ex- 
plained in detail. In section ^, we discuss the metric perturbation in the external 
scheme. In section ^, the equations of motion for a monopole particle are derived 
by using the method (1). The method (1) is applied only to the case of a monopole 
particle because of the difficulty in constructing the perturbed metric of a Kerr black 
hole. The case for a spinning particle is considered in section ^ by using the method 

(2). 

Throughout this chapter we assume that the background metric satisfies the 
vacuum Einstein equationshil. Hence in the following calculations we use the fact 
that the background Ricci tensor vanishes; 

R^^u = 0. (M) 

§2. Matched Asymptotic Expansion 

The matched asymptotic expansion is a technique with which the same physical 
quantities derived in different zones by two different approximation schemes are 
matched in the overlapping region to obtain an approximate solution valid in the 
whole region. We first prepare the metrics in both internal and external zones by 
using different approximation schemes. The internal zone is the region where the 
self-gravity of the particle dominates while the external zone is the region where the 
background geometry dominates the full geometry. 

In the internal zone, we assume that the metric can be described by that of a 
black hole plus perturbation. Namely, we assume that the particle is represented by 
a Schwarzschild/Kerr black hole in the lowest order of approximation. In the present 
case, the perturbation is caused by the tidal effect of the curvature of the spacetime in 
which the particle travels. As mentioned in Introduction, we call this construction 
of the metric the internal scheme. In order to make this scheme valid, the linear 
extension of the internal zone around the particle must be much smaller than the 
background curvature scale L. We introduce the coordinate {A"} = {T, A*} (a = 
0,1,2,3; i = 1,2,3) for the interiiaL scheme and |A|(:= vA*A*) is assumed to 



represent the physical distance scalaiUlil. Then the internal scheme is valid when 

|A|<CL, (2-1) 

where L is the length scale of the background curvature. 



*•* The result is not altered even if we assume that the background spacetime is vacuum just 
around the particle. 

**' In this chapter, we adopt the Minkowskian summation rule on a,b,- ■ ■, and the Kronecker 
summation rule on i,j, ■ ■ ■ over the repeated indices. 
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In the external zone, we expect that the metric is well approximated by the 
perturbation induced by a point source on a given background spacetime. We call 
this construction of the metric the external scheme. This approximation scheme is 
valid when the self-gravity of the particle is sufficiently weak, that is, 

Gm < \X\ , (2-2) 

where (Gm) is the scale of Schwarzschild radius. As the point source is placed where 
the external scheme is invalid, there is no matter source in the external zone. Thus 
the external metric is given by a vacuum solution of the Einstein equations. 

We require that the metrics obtained in both schemes be matched in the overlap- 
ping region of both zones, by considering a coordinate transformation between the 
internal and external metrics. Safely, we may assume the existence of the matching 
region as long as 

Gm<^L, (2-3) 

is satisfied. For definiteness, we set the matching radius at 

\X\ ~ (GmL)i/^ (2-4) 

in the spatial coordinates of the internal scheme, X^. Then writing down the met- 
ric in the internal scheme, we have two independent small parameters \X\/L and 
Gm/\X\ in the matching region. The power expansion with respect to these two 
small parameters allows us to consider the matching order by order. 

First we consider the expansion of the internal scheme. Recalling that the per- 
turbation in the internal zone is induced by the external curvature which has a 
characteristic length scale L, the metric can be expanded in powers of |^|/./^ as 

gabiX) = (°)i?a6(X) + h^)Hab{X) + -Li^) Hab{X) + • • • , (2-5) 

where ^^'Hah{X) is the unperturbed black hole metric. We expect that ^^>Hah{X) 
will be given by the standard linear perturbation of the black hole. Later, we find 
that ^"^^ Hab{X) can be consistently set to zero, which is in accordance with the 
notion that the spacetime curvature is of 0(1/L^). Thus the standard black hole 
perturbation theory applies up to ^'^' Hab{X). Further we expand the metric with 
respect to Gm/\X\ which is also small at the matching radius: 

(0)i7,fe(X) = r/,5 + Gm^^HabiX) + (Gm)2[°j//,,(X) + • • • , 

J^ '^ab{^) - -^(o)^afel-^ j + -^(i)-nabi-^ j H ^ (2)^ab\^) H , 

-^^ 'Hab[X) - j^^QS^Hab(X) + -j^(^i)^ab(X) H j^^^^Hab[X) H (2-b) 

Note that, from the definitions of the expansion parameters, the J^-f Hab component 
of the metric behaves as 

(m) 



(n) 



Hab ~ \Xr-''. (2-7) 
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The explicit form of the coordinate transformation from the general coordinates 
of a background metric {x^} to the coordinates of the internal scheme j-'^"} will 
be discussed in section y for the method (1) and in section g for the method (2). 
Assuming the matching can be consistently done, the full metric in the external 
scheme g^u{x) is written in terms of the internal coordinates as 

gabiX)dX^dX'' = ~g^,{x)dx^'dx'' . (2-8) 

Generally, as the external metric can be expanded by Gm/\X\, we write it as 

~9ab{X) = gabiX) + Gm(i)/i,b(X) + {Gm)\2)hab{X) + ■■■. (2-9) 

Then Gmn\hab{X) can be recognized as the linear perturbation on the background 
gab{X). Further we expand it with respect to \X\/L as 

gab{X) = f^^KbiX) + ^%hab{X) + -^[o)^a6(X) + • • • , 

Gm(i)/i„b(X) = Gmf^KbiX) + ^j;)/j,,(X) + ^g/i„,(X) + • • • , 
{Gm)\^)Kb{X) = {Gmff^^Kb{X) + i^«/j,,(X) 



As before. 



+ ^g/^..(X) + .... (2-10) 



j:;^)/^,^ ~ \xr-^. (2-11) 



For brevity, we call ^\hah or /^^ Hat the (™)-component and the matching con- 
dition for them as the (^) matching. In the matching region {\X\ ~ (GmL)^''^), 
the (^)-component is of O ( {Gm/L)^"^~^^>''^ j . The matching condition requires that 

all the corresponding terms in Eqs. (2-6) and (2T0) should be identical. Then the 
matching condition is given by equating the terms of the same power in \X\ in both 
schemes to desired accuracy. Thus the condition for the (™) matching is 



{GmTim')^ _ sr ^GmY i^')rr , n (^m)' 



n+l 



V^ V-T/zt; (m'), _ Y^ VLT/zy (m'jrr , ^ K'-""') |y|(>n-n) 

m' —n'^m — n m' —n' ^m — n \ / 

(2-12) 

§3. External Scheme 

As we assume that the gravitational radius of the particle, Gm, is small compared 
with the length scale of the background curvature, L, we approximate Sg^i, by the 
linear perturbation induced by a point-like particle, h^i,, in the whole spacetime 
region except for the vicinity of the world line of the particle. The calculation is 
performed in an analogous manner to the case of the scalar and vector perturbations 
developed by DeWitt and BrehmecP. 
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We take a Green function approach to study the linear perturbation of the metric 
generated by a point source. In order to calculate the tensor Green function in a 
background covariant manner, we begin with introducing the concept of bi-tensors. 

3.1. Bi-tensor formalism 

Bi-tensors are tensors which depend on two distinct spacetime points, say, x and 
z, so that they can have two types of indices. The simplest example is given by a 
direct product of tensors at the points x and z as 

A^"^{x,z) = B^'{x)C'^{z). (3-1) 

In what follows, we use x for a field point and z for a point on the particle trajectory, 
and assign the letters a, (3, 7, 6, e, C, ?? for the tensor indices of z and fi, u, ^, p, a 
for x. 

Basic bi-tensors used in our calculations are half the squared geodetic interval 

a{x,z), 

(y{x,z) = -g^"" {x)a-^{x, z)a-^{x, z) = -g°''^ {z)a.a{x, z)a.i3{x, z) , 

lim a(x, z)-a = lim afx, z)-a = , (3-2) 

and the geodetic parallel displacement bi-vector, 

9f,aAx^ Z)g'"^{x)<y-a{x, Z) = 0, g^,a;f3{x, z)g'^''' {z)a-^{x, z) = 0, 

lim5/(x,z)=5/. (3-3) 



a;— »z 



These are used to expand bi-tensors around the orbit of a particle. For example, we 
have 

A'^{x,z) = lim (A"(x',^)-(T;.,(x,x')A"''^'(x',z)+0(e2) ) , (3-4) 

B^^ix) = r^ix, z) (B'^iz) - a,p{x, z)B^'P{z) + 0{e^)\ , (3-5) 

for a small geodetic interval between x and z, where e = y^2\a{x, z)\. These relations 
can be verified by taking the x ^ z limit of their repeated derivatives. 

By evaluating the repeated derivatives of Eqs. (3-2) and (3-3) in the coincidence 
limit x ^ z, we obtain some useful formulas for expansion in e: 

1 . „ 

a-apix, z) = gaf3{z) - -Ra^f3 {z)a-^{x, z)a-s{x, z) + 0(e ) , (3-6) 

(^■,^lp{x,z) = -g^'^{x,z) (gapiz) + -Ra'yp6iz)cr'^^ (x, z)a'\x, z)) 

+Oie') , (3-7) 

5^^^(x,z) = -^g^^x,z)R"^ps{zy''{x,z)+0{e^), 
r'^-Ax, z) = -\r^{x, z)g,\x, z)R^p^s{z)(T-'\x, z) + 0{e^) . (3-8) 
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We also introduce the van Vleck-Morette determinant, A{x,z): 

A{x,z):=\r^{z,x)a,^f3{x,z)\, (3-9) 

which appears in the expression of the Green function later. 

3.2. Tensor Green function 

We consider the linearized Einstein equations. We introduce the trace-reversed 
metric perturbation, 

i^fMuix) = h^y{x) - -g^y{x)h{x) , (3-10) 

and set the harmonic gauge condition, 

^P^'''.,{x)=0, (3-11) 

where h{x) and 'ip{x) are the trace of h^^'^{x) and that of 'ip^'^{x), respectively, and 
the semicolon means the covariant derivative with respect to the background metric. 
In this gauge, the linearized Einstein equations become 

_l^/^'^;5^(a;) _ R''^''p{x)iP^P{x) = 8ttGT'"'{x) . (3-12) 

Thus we define the tensor Green function G^^"^(x,z) which satisfies 
Gf"''"^'^.^{x, z)+2i?^5%(x)G«^°'5(x, z) 

= -25°(^(x, z)r^^(x, z)^—^^^^ , (3-13) 



-9 

where g is the determinant of the metric g^j,i,{x). 

First we consider the elementary solution G't' {x,z) which satisfies Eq. (3-13) 
except at the o"(a;, z) ^ limit and takes the Hadamard form, 

Gr"^(x,z) = ^l'^—^^+v^^--P[x,z)\og\a{x,z)\ 

+u;'^^"^(x,z) I . (3-14) 

The bi-tensors n^'^"^(x, z), v^'^"'^{x, z) and w'^'^°'^{x, z) are regular in the o"(x, z) ^ 
limit and u^^°^^{x,z) satisfies the normalization condition, 

lim «^'^"^(x, z) = lim 25"(^(x, z)f^'^{x, z) . (3-15) 

If we put the form (3-14) into the left hand side of Eq. (3-13), the terms can be classi- 
fied into three parts. One is the terms which contain the factor 1/(T^(x, z) manifestly 
and another is the terms which contain log |o"(x, z)|. The remaining terms have no 
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singular behavior at the a{x,z) — > hmit. Since the form (3-14) is redundant, we 
can set these three sets to vanish separately: 

( 2n''^°^'«(x, z) - ^-i^uf^'^'^ix, z)\ a,^{x, z) = , (3-16) 

\ A[x,z) J 

^/.^a/3;5.^(^^ Z) + 2R''^%{X)V^P'''^{X, z) = , (3-17) 

2v^'''''^{x,z) + l2v^'''''^'^x,z) - ^-j^v^""'^{x,z)] a..{x,z) 
V A{x,z) J '^ 

+u'"''''^'^.^{x, z) + 2i?^g%(x)n«''"^(x, z) 

+ (w^"''"^'^,i:{x, z) + 2R^'i:%{x)w^P'"^{x, z)) a{x, z) = . (3-18) 

Equation (3-16) is solved with the normalization (3-15) as 



u'"''"^{x,z) = 2g''^P{x,z)f'>'^{x,z)^A{x,z). (3-19) 

The bi-tensors vP''^'^^{x, z) and wP''^"^{x, z) are to be determined by solving Eqs. (3-17) 
and (3-18). The bi-tensor w^'^°^^{x,z) is not needed but the bi-tensor v^'^°'^{x,z) 
plays an important role in the following discussion. Although it is difficult to find 
the solution of vP^°'^{x,z) in an arbitrary background spacetime, its explicit form 
is not required for the succeeding discussions. However it is important to note that 
yi^'^'^p (^x , z) is uniquely determined. The reason is as follows. Prom Eq. (3T7) one 
finds it satisfies a hyperbolic equation. Hence the problem is if its Cauchy data are 
unique or not. First we note the coincidence limit of Eq. (3-18), which gives 

h_mz;^^°^(x,z) = lim25"(5(z,x)/p)(z,x)i?'^«'^^(x). (3-20) 

Then taking the null limit a{x,z) ^ of Eq. (3-18), we obtain the first order 
differential equation for vP''^"^{x, z) which can be solved along a null geodesic. Thus 
this equation with the boundary condition (3-20) uniquely determines v^'^°^^{x, z) 
on the light cone emanating from z. Therefore the hyperbolic equation (3-17) has a 
unique solution. We also mention that v^^^°^^{x,z) is divergence free, 

v^"''^^,u(x,z) = ^. (3-21) 

To see this we note the harmonic gauge condition on the Green function requires 

limt.^^°^.^(x,z) =0. (3-22) 

We also see that the equation for v'^'^'^^-,y{x, z) follows from Eq. (3-17), 

■t;^'^"^;,(x,z)]'«;g = 0, (3-23) 

where we have used the fact R^^'^^-p = 0, which is proved by contracting the Bianchi 
identities for the vacuum case. Thus we conclude that Eq. (3-21) holds everywhere. 
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The Feynman propagator G^ {x, z) can be derived from the elementary solu- 
tion Gt^ {x,z) by the ze-prescription. 

^M^»/3. ._ 1 l u^^-P{x,z) 
^ ^ ' '~ {2^Y\a{x,z)+ie 

+7;^''"^(x, z) log(a(x, z) + ie) + w^""'l^{x, z) j . (3-24) 

The imaginary part of the Feynman propagator Gp' {x,z) gives the symmetric 
Green function G^^°^^{x,z), from which we can obtain the retarded Green function 
G'^^^ {x,z), and the advanced Green function G'^^ {x,z) as 

G^"''"^{x,z) = --Im \g'P"'^{x,z) 

= — \u^'''''^(x,z)5(a(x,z))-v^"''''^(x,z)9(-a(x,z))] , (3-25) 
Svr L J 

GZf{^,z) = 2e[E{x),z]G^'''^''{x,z), (3-26) 

GZ:^ix,z) = 2e[z,Eix)]&'^'^^ix,z), (3-27) 

where U{x) is an arbitrary space-like hypersurface containing x, and 9[U{x),z] = 
1 — 6[z, X'(x)] is equal to 1 when z lies in the past of S{x) and vanishes when z lies 
in the future. 

3.3. Metric perturbation 

Using the above obtained retarded Green function, we compute the trace-reversed 
metric perturbation ip^^[x) induced by a point-like particle. We assume the energy- 
momentum tensor of the form, 

^ — ^ (mono) "T -f (spin) ' y"^ ^°^ 

TlZono)i^) = m f dTv^{x,T)v'^{x,T)^-!-^^^^^ , (3-29) 

^g^n) = -mJdTV^ (s^('^{x,TX\x,t /'^^^^^^^A , (3-30) 

v^^{x,T)=ra{x,z{T))i^{T), (3-31) 

S^''ix,T)=raix,ziT))rpix,z{T))S"^iT), (3-32) 

where z"'{T) = dz'^/dT, m is the mass of the particle and S'^^{T) is an anti- 
symmetric tensor representing the specific spin of the particle per unit mass. We 
call it the spin tensor of the particle and assume that it satisfies the center of mass 
condition, 

5„/3(r)i^(T) = 0. (3-33) 

In Chapter lj_section 11, we have given the energy-momentum tensor of a spin- 
ning test particle .til There the four- velocity of the orbit v" = i" is distinguished from 

*' Note that Sap there corresponds to mSap here. 
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the specific four-momentuiii of the particle u" = p'^/m. The difference is 0{S'^ /L?') 
where S is the magnitude of the spin tensor S :- 



' Sa/3S"^/2. Here we ignore this 
difference because of the following reason. Since the particle is assumed to represent 
a black hole, m will be identified with the black hole mass and S with the Kerr spin 
parameter a. Therefore S is assumed to be of order Gm, hence the difference be- 
tween f" and m" is 0{{Gm/L)'^). Since we are interested in the radiation reaction of 
0{Gm/L'^) to the equations of motion, we may consistently neglect this difference. 
At this point, we must comment on the reason why we may assume the point-like 
particle for the source. Even in the linear perturbation, in order to generate a general 
gravitational field in the external zone, we need to consider a source with arbitrary 
higher multipole moments. tJ However, the £-th moment of the gravitational field will 
be 0{{Gm/\X\Y~^^) in the matching region if the particle represents a black hole. 
As we shall see in the following discussions, we find it is not necessary to consider 
the matchings at 0( (Gm)^j or higher in order to derive the equations of motion 

with the reaction force of 0{Gm/L?'). Hence the moments higher than the spin can 
be consistently neglected. 

We should also note that the metric perturbation induced by 7'/^pj„^ is OuGm)'^ j 
for S = 0{Gm). At first glance, one might think that this implies the necessity of 
the second order perturbation theory if we are to incorporate the spin effect of 
the particle in the expansion with respect to Gm in a consistent way. However, 
provided that the construction of the metric by the matched asymptotic expansion 
is consistent, the second order perturbation theory turns out to be unnecessary. In 
fact, we shall find that the spin-induced metric perturbation of Of (Gm)^ j gives rise 
to the leading order spin-curvature coupling term of 0{Gm/Lp') in the equations 
of motion, while the spin-independent metric perturbation of 0( (Gm)^) does not 

contribute to the reaction force term at 0{Gm/L'^). 

Without any further approximation, the metric perturbation due to the point- 
like particle becomes 



'^^''(x) = 2Gm\ 



1 



+ 



&{x,z{T)) 
a{x,z(T)) 






&Hx,ziT)) 

+-^^^^u^\pni^,z{T))S^'^{T)z^{T) 
1 d 



a'^{x,z{T))dT 
1 



[u>^\p{x,z{T))a.,^{x,z{T))S^^{T)i''{T) 



+- 



&{x,z{T)) 



v^^\f3{x,ziT))a,^ix,ziT))S^^{T)i''{T) 



^T=TR,tix) 



*' A distributional form o£,the energy-momentum tensor with arbitrary higher multipole mo- 
ments was discussed by DixonEr. 
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V 

+v^\fS,^{x,z{T))S''^{T)i''{T)\ J , (3-34) 

where Tfiet{x) is the retarded time of the particle and is a scalar function which is 
determined by 

a(x,z(TRei)) = 0, e{S{x),z{TRet)) = l. (3-35) 

Since the retarded time Tuetix) is not convenient for specifying the field point 
X around the particle trajectory in the following computations, we introduce a new 
specification of x as follows. We foliate the spacetime with spacelike 3-surfaces 
perpendicular to the particle trajectory. Specifically, the 3-surfaces are defined as a 
one-parameter family of T by the relation, o".„(x, z(T))i"(T) = 0. We denote the 
value of T of the 3-surface containing the point x by Te. That is 

a;„(x,z(r,))i°(T,) = 0, (3-36) 

where we have introduced the notation, 

Q-a(.X,z{T^)) := [Q;a(x,z)]^=^(T,) , 

Q-f,{x, z{T^)) := [Q-i,{x, z%=^(T^) . (3-37) 

Note that 

[g(x, z{T^))].^ = Q.,^{x, z{T^)) + Q;„(x, z(T,))i°(T,)r,;^ . (3-38) 

We use a-^a{x-,z{Tx)) to distinguish the spatial points on the same 3-surface, and 
denote the spatial distance from z{Tx) to x by 

e{x):=^2a{x,z{T^)). (3-39) 

In the matching region, we have 

Gm < e{x) < L. (3-40) 

To obtain the external metric in the matching region, we first consider the e- 
expansion of the time retardation, 5iiet{x), 

5Ret{x) := TRet{x) - T, . (3-41) 

It is given by expanding Eq. (3-36) as 

0= [a(x,z(r))]^^y^^^(^) 
= a{x, z{Tx)) + o-(x, z{Tx))5Ret{x) 

+-a{x,z{T^))5\^t{x) + ^ o^ {x^z{T^))6\^t{x) 
+l-^-(x,z(r,))4,,(x)+0(e5). (3-42) 
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Using Eqs. (3-6), (3-36), and the normalization condition, (dz/dT)'^ = —l+0{Gni/L), 
which will be proved to be consistent later, each term in the above is computed as 



:e^(x) 



a{x,z(Tx)) 

a{x,z{T^)) = a.a{x,z{T^))z"'{T^) =0, 
a{x,z(Tx)) =: -k^{x) 

= a;,^(x,z(T,))i"(T,)i'3(r,) + (j;„(x,z(T,))r(T,) 



(3-43) 
(3-44) 



9aMTx)) - ^Ra''pH<T^))^.^{^,4T,))a.5{x,z{T,))] i"(T,)i^(r,) 



+a;„(a;,z(r,))r(T,) + 0(e3), 
a {x, z{T^)) = a.a{x, z{T^)) 'z "(T,) + 0{e^) , 
a {x,z{T^)) = -ga(s{z{Tx)yz"{T^)z''{T^) + 0{e) . 



(3-45) 
(3-46) 
(3-47) 



where we have introduced k{x) to denote y/—a{x, z{Tx))- From these, we obtain 



5Ret{x) = -e{x)K'\x){l - -e{x)n-^{x) -z "{T^)a.a{x,z{T^)) 



24 



e\x)^-\x)zHT,)]+Oie^). 



(3-48) 



With the help of Eq. (3-48), we then obtain the expansion of various terms in 
Eq. (3-34). We have 



1 



&{x,z{T)) 



T=TB.et{x) 



= T^ ( 1 - \<x) -z "(r.)a;«(x, z{T,)) - Je2(x)z2(r.) + oie'^m) 
e[x)K{x) V 3 8 / 

In order to obtain the expansion of u^^'^°'"{x, z) given by Eq. (3-19), we also need the 
following expansions: 



A^/\x,z{T)) 



1 + 0(e3 



(3-50) 



T=TR^t(x) 



g^a{x,z{T)) 



T=TR,t{x) 



g^a{x,z{T^)) - -g/ {x, z{T^))Rc,p^5{z{T^))a'^ {x, z{T^))z^ {T^)e{x) 



i"(r) 



+0{e') , 



T=Tii,t(x) 



(3-51) 
z'^iT^) - e{x)K-\x)z"{T,) + ^e\x) z "(T,) + 0{^-52) 



In the above expressions there appear higher derivatives of i, such as z and z, 
where a dot means the covariant derivative D/dT along the trajectory of the particle. 
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Since we are considering the case in which the radiation reaction force is 0{Gm/Lp'), 
it is reasonable to assume these derivatives are smaher by a factor of 0{1/Tr)^ i.e., 

"";'-f)^_L_<o(iW) („>i), (3.53) 

dT^+i TrU'-^ V^"+V ^ - ^' ^ ^ 

where T^ = 0{L'^ /{Gm)) is the reaction time scale. We shall find that this is 
consistent with the equations of motion in the end. 

Keeping this fact in mind, and using Eqs. (3-49) ~ (3-52), we obtain the e- 
expansion of the trace-reversed metric perturbation, Eq. (3-34), as 

r'' = V'fi„,) + ^f;,„) + ^sio , (3-54) 

where 



_4^{a^/3) + 2z<a^^z'R^&,^'^z^^ - 2eR''/si^i^ + 0{e^) \ (3-55) 



C«„^(^) = -"^Gm-g^r, -^i("5^)V;, + Oi{Gm)e') , (3-56) 



X I -£^drYt;"%,(z(r,),z(r'))i°'(rOi^'(r') 



"/3 ,., -(^(T^),z{T'))S^''''{T')if^'{T') 



+^"^'/3';y(^(T.),z(r'))5^" (r')i^ (TO 



J —oo \ 

+^"^',3';y^^(^(r.),z(T'))5^'°'(T')i^'(T0')+O(e2()3|57) 

where ^^ = g^(x, ^(Tj:)). The part tpu'^^) is called the tail term because it is not due 
to the direct light cone propagation of waves but due to multiple curvature scattering 
of waves as described by the v^^^'^^{x, z) term in the Green function. 

3.4. Transformation to the internal coordinates 

In order to write down the external metric in terms of the internal coordinates, 
we consider a coordinate transformation from x to {^"} given in the form, 

a.,c,{x,z{T)) = -F^{T,X). (3-58) 

We restrict our consideration on a coordinate transformation which satisfies the 
following requirements. We assume X^ = corresponds to the center of the particle, 
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x" = z°^{T), hence J^a = at X* = 0. We also assume that the right hand side of 
Eq. (3-58) can be expanded in positive powers of X^ as 

J^a{T,X) = fai{T)X' + ^Uij{T)X'X^ + ^f^ijk{T)X'X^X'' + ■■■ .(3-59) 

Although it is possible that there appear more complicated terms such as X^X^ /\X\, 
we simply ignore such kinds of terms. We shall find it is consistent within the order 
of the approximation to which we are going to develop our consideration. Here 
faix--in{T) is totally symmetric for ii • • • i„ and is at most of 0{L~'^'^~^>). Using 
Eqs. (3-6) and (3-7), the total derivative of Eq. (3-58) gives the important relation, 

dz°' Df°'- 1 D f " • • 

rMT),x)dx'^ = I ^(T) + -j^{T)X^ + -^{T)X^X^ 



lR%,s{ziT))f,{T)^{T)f',{T)X'xn 



dT 



dT 






k 



--^R"p,s{z{T))f,{T)f\{T)fk{T)X^X' | dX^ 

+0{\Xf) (3-60) 

In the following sections, we write down the external metric in terms of the 
internal coordinates in the matching region to obtain the equations of motion. 

§4. Equations of motion for a monopole particle 

In this section, we adopt the method (1) mentioned in section |l| to derive the 
equations of motion. We restrict our consideration to the case of a monopole parti- 
cle, which is necessary because we use a well-established method to decompose the 
metric in the internal scheme by the tensor harmonics. The tensor harmonics are 
classified by the total angular momentum, J, reflecting the spherical symmetry of 
the Schwarzschild black hole. 

In the internal scheme, the monopole mode (J = 0) corresponds to the mass 
perturbation. Thus we may set this mode to zero since it is natural to suppose that 
the change of mass due to the radiation reaction is negligible. The dipole modes 
(J = 1) are related to the translation and rotation. The translation modes are 
purely gauge and thus we set them to zero to flx the center of the black hole. As 
we are considering a non-rotating black hole, we also set the notational modes to 
zero. In general, the higher modes contain gauge degrees of freedom as well as the 
physical ones. However, for these higher modes, we do not give any principle to fix 
the gauge for the moment. 
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Before the explicit computation of the (™) matching condition, we briefly review 
the construction of the scalar and vector harmonics in terms of the symmetric trace- 
free (STF) tensor!) . 

4.1. Spherical harmonics expansion 
We introduce the notation, 

^<iii2-ie>j (4'1) 

to represent the totally symmetric trace- free part of Ai^i^...i^. More explicitly in the 
cases of £ = 2, 3, 



*-<ijk> 



■^{ijk) ~ T \^ij^{km.m) "T ^jk-'^{imm) + "A;i^(jfmm) ) • l^'^j 



The spherical harmonics expansion of a scalar function A on the unit-sphere can 
be written as 

oo 

A = Y^ A<,,i,...i,>n<''n'^ ■ ■ ■ n'^> , (4-3) 

1=0 

where n* = X*/|X|. In this case, the order i, which is associated with the angular 
dependence, is equivalent to the total angular momentum, J. Thus the J mode of 
the (TT)- component of the metric perturbation is totally determined by its angular 
dependence. Namely, the terms in the (TT)-component of the metric perturbation 
which contain 

1, n\ n<'n^>, (4-4) 

correspond to the J = 0, 1, 2 modes, respectively. 

Next we consider the expansion of a vector field Ai, 

oo 

Ai = Y, Ai<,,i,...,,>n<''v}^ ■ ■ ■ n'^>. (4-5) 

e=o 

In this case the term of the i-th order in the angular dependence is decomposed 
into J = i + 1, i and i — 1. This is done by using the Clebsch-Gordan reduction 
formulaEI), 

TIT- - »(+) ^—^^ rW I ^^"^ ^ o(-) (^f^^ 

where Ti-^i^.-.i^ is a STF tensor of order i and 



11Z2---11+1 
o(0) 
ili2---ii 
o(-) 



'-^ <il+l^iii2---ii>J 

'-^ j -^ k<iii2---ii—i^il>jkj 

U j''-jiii2---ii-\- (4''j 



We perform the decomposition explicitly for H. < 2 here. For £ = 0, there exists 
no J = mode and it trivially corresponds to the J = 1 mode. For i = 1, the 
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decomposition is performed as 

1 \ 11. 

(4-8) 



AijTi 



^fe') ~ '^^ij^kk j + A[ij] + l^^ijAkk 



and the first, second and third terms in the square brackets correspond to the J = 2 
,1 and modes, respectively. For £ = 2, we obtain the decomposition formula as 



where 



A ■; + ^f ■ -B^^) +^r^- •R(^) 



n<%'^>, (4-9) 



?{2) 



1 



-Ojj r) V^fc<m«>^ifem I ^fc<mjr'>^i 



tj 2 



ikm) 



bI^'^ = Ai<Jk>S^, , (4-10) 

and the first, second and third terms correspond to the J = 3, 2 and 1 modes, 
respectively. 

We omit the general discussion on the expansion of the tensor field and we shall 
give a specific argument when necessary. 

4.2. Geodesies; (q) and (g) matehing 

We begin with the (g) and (J) matchings which are, respectively, of 0{{Gm/L)^) 
and of 0{{Gm/LY''^) in the matching region. First we consider the external scheme. 
In these matchings the external metric is the background itself. Here, the necessary 
order of expansion in \X\ is 0(|X|). We note 

g^^{x)dx^dx'' = gc,p{z)g''^{z, x)f^{z, x)dx^dx'' . (4-11) 

Then from Eq. (3-60), we get 

g,Ax)dx^dx^ = ((^)' (T) + 2^iT)^{T)X^^ dT^ 

{ dz'^ dz"" ■ Df A 

+2 [^{T)fUT) + —{T)U,{T)X^ + f-,{T)^{T)xA dTdX^ 



+ [nmU.iT) + 2r,{T)f^,k{T)X') dX^dX^ 

+0 [m . (4.12) 



Comparing the above equation with Eq. (2-10) and looking at the dependence on X, 

(0), ^. (1), 



one can readily extract out (nlhab and /nlhab to the lowest order in Gm/L. 



Next we consider the internal scheme. The ([D-component is trivially given by 
the flat Minkowski metric. Hence the (q) matching becomes 

"^ = (^)' ^^^ + ^ (^) ' (^^)-component, (4-13) 

= ^{T)U,{T) + O (^^j , (Ti)-component, (4-14) 

S^, = r^{T)faJ{T) + O (^) , (ii)-component. (4-15) 
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Equations (4-14) and (4-15) indicate that /"*(T) are spatial triad basis along the 
orbit, i.e., 

rk{T)fMT) = g-^z{T)) + ^m^iT) + O (^) . (4-16) 

To know the (Q)-component of the internal scheme, it is better to consider all 
the (^)-coinponents at the same time. Namely we consider the linear perturbation 
of the black hole ^^> Hab- For this purpose, we consider the harmonic decomposition 
of linear perturbation as discussed in subsection 4.1. Since the time scale associated 
with the perturbation should be of the order of the background curvature scale L, it 
is much larger than the matching radius {GmLY''^ . Therefore the perturbation may 
be regarded as static. It is known that all the physical static perturbations regular 
on the black hole horizon behave as ~ \X\ asymptotically where J is the angular 
momentum eigenvalue. However, in ,{Hab, there exists no term which behaves as 

~ \X\"^ {m > 2). Hence, except for gauge degrees of freedom, SJ^Hab contain only 
J = 0, 1 modes. As mentioned before, we set the perturbation of J = 0, 1 modes to 
zero. Thus we conclude that we may set 

Si5^ab = 0, (4-17) 

for all n. This is the gauge condition we adopt for the internal scheme at 0{1/L). 
In particular this condition gives the (q) matching as 

° = ^^(^)^(^)^^ + ^ (^?) ' iTT)-co^voneni, (4-18) 






' Gm \X\ 
+0 ( — — ^ ) , (Ti)-component, (4-19) 



= 2f^u{T)rj)k{T)X^ + O -—'-^ , (ii)-component. (4-20) 



^k ^ f Gm \X\ 
Ta{i\J^ )J j)ky-L )^ "^ l~L — L 

Then the covariant T-derivative of Eq. (4T3) and that of Eq. (4-14) with Eq. (4T8) 
result in the background geodetic motion. 

One can see from Eq. (4-13) that the internal time coordinate T becomes a proper 
time of the orbit in the lowest order in Gra/L. In the same manner, Eq. (4-18) and 
the covariant T-derivative of Eq. (4T5) with (ij)-antisymmetric part of Eq. (4T9) 
give the geodetic parallel transport of the triad /"^(T), 

Further, from Eqs. (4-19) and (4-20), we can see 

/%(r) = o(^i). (4-23) 
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4.3. Hypersurface condition; (g) matching 

We now proceed to the (g) matching, in which the external metric is stih given 
by the background but there appear non-trivial perturbations in the internal scheme. 
Although it is of 0{Gm/L) in the matching region and 0{{Gm/LY''^) higher than 
the remaining (5)-component, we consider it first for the reason which will be clarified 
below. 

(2) 

In order to obtain )Q\hab, we expand the external metric in terms of the internal 
coordinates up to 0(|Xp), i.e., we have to go one order higher than Eq. (4-12). Then 
the (g) matching becomes 



(2) 



dz^ ,„, „fl ,„, dz'^ 



j^\;\Htt = -Ra(^,5{z{T)) — {T)fMT) — {T)f,{T)X^X^ 

— — J-^- , (rr)-component, (4-24) 

2 d7^ 

+0 —j-^-jj- , (ri)-component, (4-25) 

^(o)-^«i ^ faii{T)f°'j)ki{T)X X 

-^-R^p,s{z{T))r,{T)f^k{T)P,{T)fi{T)X^X' 

^(Gm\X\'^\ , ^ , ^ 

+0 — —T- , (ijj-component, (4-26) 



where Eqs. (4-22) and (4-23) have been used to simplify the expressions. Since we 
have set ,{Hah = 0, the first non-trivial perturbations of the internal metric appear 

(2) 

in ){Hab- Hence they describe the linear perturbation of the black hole metric in 
the internal scheme. Then we have to fix the gauge condition for this perturbation 

(2) 

to perform the matching. For )J.Hab, since the physical perturbation contained in 
it is quadrupolar, we fix the gauge so that all the J modes except J = 2 are zero. 
Then the (g) matching becomes as follows. 

First consider the (rT)-component of the metric. The right hand side of Eq. (4-24) 
may contain J = 0, 2 modes. The J = mode, however, vanishes because of the 
background Ricci flatness. Hence this matching just determines the physical pertur- 
bation in the (TT)-component. 

As for the (ri)-component, the right hand side of Eq. (4-25) may contain J = 1, 
2, 3 modes. As before, the J = 2 mode just determines the physical perturbation of 
the (Ti)-component. So we put J = 0, 3 modes to zero. However, they are found 
to be absent in the second term of Eq. (4-25). To see this we first decompose its 
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angular dependence, 



-^Ra/S^Spi {f<jf k>X^^X > + -f^kf k\X\'^ 



(4-27) 



Using Eq. ( 4-16 ) and the fact that the Ricci tensor vanishes, the second term in the 
parentheses is rewritten as 






3 dT 



dT dT' 



(4-28) 



and is found to be zero due to the symmetry of the Riemann tensor. The first term 
in the parentheses of Eq. (4-27) is decomposed further with the aid of the formulas 
(4-10) and {^ as 



dz^ 



n .fP.f^, 



(2)7/35 






+ 7^«<i^fc> 



,p{ihl3s\ x<^X^> 



(4-29) 



where 



p(2)7/35 \ ( fl fP f^ c -L- n fP f^ p 

'^ ij • — ry \i Taj <nJ iy'^jmn ~r J mj <nj jy'^irr 

fI'^'^' := \ [Puf^fk + PkPkf\) - \p^f^fk . 



(4-30) 



It is easy to see that the first and third terms in the parentheses of Eq. ( 4-291 ) vanish 
due to the symmetry of the Riemann tensor and the Ricci flatness. Thus only the 
J = 2 mode remains in the second term in the right hand side of Eq. (4-25). 

Decomposing the first term in the right hand side of Eq. (4-25) in a similar 
manner, we find it contains J = 1, 3 modes as well as J = 2 mode. Putting the 
J = 1 mode to zero, we obtain 



^(T)/,<„(r) = o i—-. 



Putting the J = 3 mode to zero gives 
1 dz 



2 dT 



{T)f^<,,k>{T)XiX^ = 



'Gm\xf 
L L2 



(4-31) 



(4-32) 



Then combining this with Eq. ( 4-31| ), we find 



dz"' ,„, „ ,^, ^ f Gm 1 



From Eqs. (4-14), (4-23) and (^331), we find 



(4-33) 



dz°' dz"" 

^{T)a,^ ix{T,X),z{T)) = - — {T)J^^{T,X) = O 



\X\' 
1? 



-L 



(4-34) 



to the lowest order in Gm/L. Comparing this with the hypersurface condition 
of Tx-, Eq. (3-36), one finds that the T = constant hypersurface differs from the 
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Tx = constant hypersurface only by O(e^) = 0(|X|^). It then follows that all the 
calculations done in section y remain valid even if we replace Eq. (3-36) with 

a.aix, z(r,))i°(r,) = O(eVi'). (4-35) 

Thus T can be identified Tx to the lowest order in Gm/L. The reason why we 
have done the (q) matching prior to the remaining (5) matching is to establish this 
equivalence of T and T^. 

Turning to the (ij)-component, it may contain J = ~ 4 modes, we first note 
that the second term of Eq. (4-26) contains only J = 2 mode. This can be seen as 
follows. First, we define the spatial triad components of the Riemann tensor by 

Rijkm '■= RaP^sf^if jPkf m ■ (4-36) 

Introducing a symmetric tensor defined by 
we can express R^'^'J'^ in terms of TZij as 

^ijkm 

Then the symmetric tensor TZij is decomposed into STF tensors as 



Ri.km = e™^■e^'^"^7^„, . (4-38) 



TZij = T^<ij> + i^^ijT^kk- (4-39) 

Counting the number of indices, we find that the first and second terms in Eq. ( [4 -391 ) 
correspond to J = 2 and modes, respectively. However, again owing to the sym- 
metry of the Riemann tensor and the Ricci flatness, the J = mode vanishes and 
only the J = 2 mode remains. Therefore the gauge condition for the (ij)-component 
implies 

LdT)r,),i{T)\ ^,=0 (^^) , (4-40) 



J^2 \ L L'^ 



where [• • •] j^2 means the J / 2 parts of the quantity. This will be used in the (f) 
matching below. 

4.4. External perturbation; (5) matching 

Now we proceed to the first non-trivial order in Gm/\X\. For this purpose, we 
must develop the external scheme. However, since the time slicing by the internal 
time coordinate T is now identical to that by T^ in the lowest order in Gm/L, we 
can use the previously obtained formula ( 3-541 ) for the external metric perturbation. 



Among the matchings which becomes of 0{{Gm/L)^'^) in the matching region, 
there remains the (5) matching. This matching relates the masses of the particle in 
both schemes. Since this matching is independent of L, we may regard the back- 
ground external metric as if it were flat. As is well-known, the linear perturbation 
induced by a point-like particle of mass m in the flat background spacetime is exactly 
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equal to the asymptotic metric of a Schwarzschild black hole of mass m in the linear 
order in m. This fact indicates that the matching gives a consistency condition at 
this order. 

In order to directly check the consistency, we rewrite Eq. (3-54) in terms of the 



internal coordinates. Since ,-,{h„h ~ \X\ ^, we have only to consider the -0^^ ^ term 

(1) uo I I ; J T^ [mono) 



of Eq. ( |3^ . Using Eqs. Q, (3-60) and the fact that e = y/FjT\X)F^(J\X) , 
we find 



Gm\'^>hab{X)dX''dX'> = Gm ^-AT^ + —-dX'dX' , (441) 

which corresponds to the asymptotic form of the Schwarzschild black hole of mass 
m in the harmonic coordinates. 

4.5. Radiation reaction; (}) and (\) matchings 

There are many components which become of 0{Gm/L) and 0{{Gm/L)^''^) in 
the matching region. However, we are interested in the leading order correction to the 
equations of motion with respect to Gm/L and we found in the ([]) and (g) matchings 
that in the lowest order the terms which behave as ~ IJ^I*^ or \X\^ determines the 
motion of the particle. Therefore we consider the (}) and (f ) matchings here. 

In order to perform the (J) and (f) matchings, the calculation we have done to 
obtain Eq. (4-41) must be extended to the linear order in \X\. Then the (}) matching 
equations are found as 

+0 i l-f^j ] , (TT)-component, (4-42) 

+0 i ( -p ) ) , (ri)-component, (4-43) 



L (1) 



+Ol [— ^) J , (ii)-component, (4-44) 



and the (f ) matching as 

-jYii)HTT-2—{T)^^{T)X 



~dT^ '~dT 

-^fai{T)rMT)X'X^X'X 



+Gm{ —{T) — {T)f\{T)0^p,{T)X' 
/_^ „n, /_< ,,i ,,^' ,,fc -^i 
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+^^Ra(S,5{z{T))^{T)fMT)^{T)f,{T)X'X^ 
+ol(-p-j Ml' (Tr)-component, (4-45) 

— (r)/f(r)/7(T)0„;3,(T)x^- 

+2R^p,s{z{T))^{T)f,iT)^{T)f',iT)X^ 
+^^Rap,s{z{T)) — {T)ff^,{T)f\{T)f\{T)X^X' 

+0[(-p-) Ml' (ri)-component, (4-46) 



where 



with 



Gm0afs{T) := h,^taii)a(^{z{T)), 
Gm0ap^{T) := \tail)al3;^{z{T)) , (4-47) 

h{tail)piv{x) ■■= i'(tail)f,uix) - ^9iJ,u{x)'ll^(taU){x)- (448) 



Note that /i(iaiOM'^(^) ^^ ^^^ metric perturbation due to v^^ai3{x,z) in the Green 
function. The (ij)-component of the (f ) matching is not presented here since it will 
not be used in the following discussion. 

As we have discussed in subsection 4.2, we require nlHab = 0. Thus the right 

(2) 

hand sides of Eqs. (4-42), (4-43) and (4-44) must vanish. As for )-^iHab, following the 
discussion in subsection 4.3, we set all the modes except J = 2 to zero. Inspection 
of the right hand sides of Eqs. (4-45) and (4-46) reveals that the terms involving 
the Riemann tensor are in the same forms as those appeared in Eqs. (4-24), (4-25) 
and (4-26). Hence they contain only J = 2 modes and do not give any matching 
condition. Furthermore, from Eq. (|4-40D , all the modes except J = 2 contained in 
the term involving f^jki in Eq. (4-45) vanish at the lowest order in Gm/L. Hence 
we only have to consider the remaining terms in Eqs. (4-45) and (4-46). The J = 1 
modes are extracted out to give 

dz" Df ■ dz" dzf^ 

= 2^m^m + Gm ^(r)^(r)r.(r)0„^,(r) 

+0 { — — ) y , (Tr)-component, (4-49) 
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■^(T)+Gm0„,,(r)^^'rw/3. 



= /,[,(r)^-^(T) + Gm0,^^(r)— (T)/^[,(r)r,](r) 



+oi(-P^j -j, (Ti)-component. (4-50) 

The J = mode is absent in the (TT)-coniponent, while that in the (Ti)-component 
exists but it just gives the equation which determines {dz°^ / dT) fau to the first order 
in Gm/L. 

Taking the covariant T-derivative of Eqs. (4-42) and (4-43) and using Eq. (4-49), 
we obtain the equations of motion with the 0{Gm/L'^) correction due to the radiation 
reaction, 

~iT) = -^ (e",,(T) + e-.„iT) - o^'iT)) ^(r)f (D 

Similarly the 0{Gm/L^) correction to the evolution equations of the 'triad' basis, 
f"i(T), are obtained from the covariant T-derivative of Eq. (4-44), and Eqs. (4-49) 
and (4-50). The result is 

— r,(r) = -— (0Vm + &%^iT) - Op,%T)) f,{T) — {T) 

Since the internal time coordinate T is not properly normalized in the external 
metric, we define the proper time, r = t(T), such that (dz/dr)'^ = —1. It is easy to 
see that we should choose 

^ = 1 + — e<..(T) — (T) — (T) + 0^(— ) j . (4.53) 

Since the second term on the right hand side of this equation is proportional to the 
small perturbation induced by the particle, it is guaranteed to stay small even after 
a long time interval compared with the reaction time scale Tj. = O ({Gm/L)^^L). 
Then Eq. (4-51) becomes 

D dz" ^ , 

^ _Gm f dz^dz^d^dz^ ^ ^ ^p dz^dz^ _ a5/^^^^.^0 r. 

2 \ dr dr dr dr dr dr dr dr I "^ 

.o((^)^i). 

Also, the triad basis are not properly normalized in the external metric. Thus we 
define e"j(r) as 

eai(r)e"j(r) = 5ij , (4-55) 
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e". 



:(r) = {Sij + Sij)rj - Gm{dz''/dT){dzl'/dT)p,ep^ , (4-56) 



where Sij is of 0{Gm/L) and recalling Eq. (4-43) the last term is added so as to 
satisfy the orthonormal condition, 

eai(T)((iz"/dT)(r) = 0. (4-57) 

From Eq. (4-55) we find 



s 



«j 



Gm^ _,,^.a.. w/3v,^ , ^ U '^"^^^ 



-e^p{r)f\{T)f,{T) + O ( — j . (4-58) 



Again the correction terms in ef are guaranteed to stay small. Then the evolution 
equations of the normalized triad e°'i{T) become 

§5. Equations of motion for a spinning particle 

In this section, we consider the equations of motion for a spinning particle. Dif- 
ferent from the Schwarzschild case, we cannot make use of the mode decomposition 
by the spherical harmonics since the background in the internal scheme does not 
have the spherical symmetry. Therefore, it is quite unclear for us how to fix the 
gauge in the internal scheme, and hence we cannot derive the equations of motion 
by the consistency condition of matching. 

Instead, we here apply the laws of motion and precession discussed by Thorne 
and HartleLF. As noted in section |l|, assuming the consistency between the internal 
and external schemes, we can make use of the matching condition to obtain the 
internal metric from the knowledge of the external metric. The problem to derive 
the equations of motion for a spinning particle was discussed by Thorne and HartleQ) 
and the spin-induced force was derived. The discussion given below is an extension 
of Ref. 0) in the sense that we take into account the effect of radiation reaction to 
the motion. Both derivations of the radiation reaction and the spin-induced force 
are discussed in a unified manner. 

5.1. Laws of motion and precession 

The laws of motion and precession 0) are derived from the integral identities 
given in terms of the Landau-Lifshitz pseudo-tensor, i^-L' ^^*^ ^^^ Landau-Lifshitz 
super-potential, H2_i ■ The Einstein equations can be put into the form. 



HrJ'L,,s = WnGi-g) (r-^ + tf J , (5-1) 
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where 

^a7/35^ga/3g75_ga5g7/3^ (5-2) 

+i (2g"^/'^ - 5-/^3^^) (25.C5'" - /'^/') g^^7g^%| > (5-3) 

g«/3 .^ (_g)i/2g-/3, (54) 

and a comma denotes the ordinary derivative. By construction, the following con- 
servation laws are satisfied: 



((-5)(r-^ + tf_J) =0. (5-5) 



,13 

Suppose that the internal metric around a Kerr black hole is calculated for a 
given trajectory of the particle. In terms of the internal metric we define 

IdTtG J\X\=r ^ 

+Hl% - Hit'l) , (5-7) 

where dPSj is the surface element of a two-sphere at |X| = r. Then by using the 
Einstein equations (5-1), we have the following integral identities: 

^P^{T,r)=f <fS,{-g)q^^{X), (5-8) 

Ul j\X\=r 

^J''{T,r) = J^ ^^<fS, (x\-g)t^;:_^{X) - X^{-gyt^{X)) . (5-9) 

These are called the laws of motion and precession. By explicitly evaluating the 
right hand sides of Eqs. (5-6), (5-7), (5-8) and (5-9), and eliminating P°'{T,r) and 
J''^{T,r) from the resulting equations, one obtains the equations of motion. 

5.2. Use of the matched asymptotic expansion 

In the present method, we construct the external metric and use the matching 
conditions to obtain the necessary components of the internal metric. The (^)- 
components of the internal metric are assumed to be given by the metric of a Kerr 
black hole. Since we do not construct the internal metric independently, there exists 
no a priori requirement for thus obtained internal metric to satisfy some specific 
gauge condition. Hence the transformation from the external coordinates to the 
internal ones can be rather arbitrarily chosen. Here, we make use of the knowledge 
we have obtained in section HI and we choose the coordinate conditions as follows. 
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We assume that the external metric is generated by the point-hke source, Eq. (3-28), 
and calculate the external metric in the matching region as in the previous section. 
In order to do so, the hypersurfaces of T = constant and Tx = constant should be 
identical to each other as given by Eq. ( [4-35| ). To satisfy this requirement, we adopt 
the coordinate transformation from j; to X in the form, 

a-aix, z{T)) + U{T)X' = 0{{Gmf/L). (5-10) 

This is satisfied by setting 

i/%- = L^rijk = 0{Gm/L) , (5-11) 

in Eq. (3-59). Note that, in the case of a monopole particle discussed in section 
^, the conditions that are required to guarantee Eq. ( 4-35D are obtained from the 



(Q)-matchings (n = 0,1,2). On the contrary, here we impose the conditions ( ^Tl 
by hand to guarantee Eq. ( [4-35| ). 

Furthermore, to determine the internal metric from the matching conditions, we 
set the (j'j)-components of the internal metric to zero: 

[^•^Hab = ^ (n = 0,l,2,---). (5-12) 

In the case of a monopole particle, we have found we can impose these conditions. 
However, in the present case, since we have imposed the coordinate condition ( |5T[l| ) 
by hand, it is not clear if a similar argument can be made to justify these conditions. 
Nevertheless, at least for n = 0, 1, we should be able to require the conditions 
( 5-12| ). This is because the spin of the black hole appears at OUGmY] or higher in 



the internal metric, hence the discussion we gave in the case of a monopole particle 
should be equally applicable to the (g) and (})-components of the metric. In fact, we 
see below that the conditions ( ^T2| ) for n = 0, 1 consistently determine the internal 
metric in the local rest frame by matching. 

First consider the background metric in the internal scheme. For convenience 
we define the trace-reversed (^)-components of the metric with respect to the flat 
Minkowski space: 

{m)fT _ (m) TT 1„ „cd{m) TT /c; i q\ 

(„) ^ab - („) ^ab - ^VabV („) H^d (5-13) 

Expanding the Kerr metric with respect to Gm, the (^)-components of the metric 
in the harmonic coordinates are found as 

(5-14) 

(5-15) 

(5-16) 



(O)o- 


= Vab, 


Gm'I^HTT 


AGm 


Gmf^Hn 


= 0, 


Gmf},H,, 


= 0, 
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iGmf>TT=^^, (5-18) 

(Gm)2g]Fr, = ^5.,X^ (5-19) 

{Gmf'S^H., = i^ (-2^., + ™) , (5.20) 

where Sij is the specific spin tensor in the local rest frame of the black hole. Then 
calculating the (q) and (]')-components of the external metric in the matching region, 
we find they are consistent with Eqs. (5-14) ~ (5-17) provided that z{T) and /"{T) 
satisfy the lowest order orthonormal conditions, Eqs. (4-13), (4-14) and (4-15). Fur- 
ther, the spin contribution to the metric, Eq. (5-19), can be reproduced from the 
external metric with the source (3-28) by the identification, 

Sai3{T) = S,,aT)fl + 0{{Gmf/L). (5-21) 

This fact indicates the consistency of using the point-particle energy momentum 
tensor (3-28) in the perturbation analysis. 

Keeping in mind the imposed conditions (5-11) and ( 5-12| ), the calculation of the 
(Q)-components of the external metric in the matching region gives 

^(T)^iT) = 0{Gm/L'), f^{T)^(T) = 0{Gm/L'). (5-22) 

As before, these equations imply that /"{T) is parallel transported along the particle 
trajectory at the lowest order. Similarly, the calculation of the (})-components of the 
external metric gives the same conditions as we have found in the previous section 
(see Eqs. (4-42) ~ (4-44)): 

Gm 
Y 

+0{{Gmf/L^), (5-23) 

z^mfaiiT) = -GmOap{T)z''{T)fi{T) 

+0{{Gmf/L^), (5-24) 

+0{{Gmf/L^), (5-25) 

where we have introduced 

0af3 = &al3 " l/2ga/39^5 = ^;^i'{tail)al3 ■ (5-26) 

These equations may be viewed as a coordinate condition on the internal time T. 
Clearly there is no inconsistency in them. 

Computation of the rest of components of the internal metric which are needed 
to derive the equations of motion is straightforward. The results are 

1 (2).-. 2 



z\T) = -1 - ^e^p{T) (g^^{z{T)) + 2i-(r)i/3(T) 



gj^TT = —R^p,s{z{T))i''{T)X^{T)P{T)X'{T) 



L2{or^ii- 3 
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+0(Gm\X\'^/L^), (5-27) 

-IgFT. = -lRap,s{z{T))n{T)X^{T)P{T)X\T) 

+0{Gm\X\yL^), (5-28) 

-1 g|^,, = -lR^fs,s{z{T))X^{T)X\T) (A(T)r,(r) + <5.,i°(T)i''(r)) 

+0{Gm\X\'^/L^), (5-29) 

-^Ro.p,s{z{T))z'^{T)X^{T)i^T)X'{T) 
+Gm0„/3^(T)i"(T)i^(T)XT(r) + 0{{Gm)^\X\/L^) , (5-30) 

+Gm0,^^(T)i"(T)/'3,(T)X^(r) + 0{{Gmf\X\/L^) , (5-31) 



-2Gmi?„^^5(^(T))(^^x^(r)x^(T) (r,(r)r,(T) - 5,,i°(r)i^(r)) 

-2|x|i"(T)/'^,(T)i^(r)/^,(r)' 

-2Gm0^^^iT)r,iT)fj{T)X^'iT) + 0((Gm) Vl/^') , (5-32) 

where X°'{T) = /"(T)X* and we have defined 

0al3-y '■= 0a/37 — 'i^/'^9al30 Sj ■ (5-33) 

5.3. Equations of motion 

Before evaluating Eqs. (5-6) and (5-8), let us first consider the equations for the 
spin, Eqs. (5-7) and (5-9). Equation (5-7) has a dimension of (mass) x (length) and 
we extract out the terms of 0{Gm?). Power counting of X shows that there will 
be contributions linear in the (2)-components of the metric and those from bilinear 
combinations of the (5)- and (5)-components of the metric. Then we obtain 

J'\T,r) = mSap{T)f'^'f' + 0{G^m^/L) + (r-dependent terms). (5-34) 

Equation (5-9) has a dimension of (massy and we extract out the terms of 0(Gm^/L) 
in the same way. Power counting of X shows that there will be contributions from 
bilinear combinations of the (5)- and (j'j-components of the metric, and we find that 
the right hand side of Eq. (5-9) vanishes: 

-—J'HT, r) = 0(G^m^/L^) + (r-dependent terms). (5-35) 

dT 
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Since the spatial triad are geodetic parallel transported in the background geometry 
to the leading order, Eqs. (5-34) and (5-35) result in 

^S'^'m = O (^) . (5-36) 

Thus in the test particle limit m ^ the spin tensor is parallel transported along 
the particle trajectory in the background geometry. 

We next consider Eqs. (5-6) and (5-8). Equation (5-6) has a dimension of (mass)^ 
and we extract out the terms of 0{m) and 0{Gm? /L). We find that there will be lin- 
ear contributions from (5)-, (2)-components of the metric, and bilinear contributions 
from pairs of (i)— and (5) -components of the metric. We obtain 

P°(r, r)=m + 0{G'^m^/L^) + (r-dependent terms), (5-37) 

P\T, r) = 0{G'^m^/L^) + (r-dependent terms). (5-38) 

Eq. (5-8) has a dimension of {mass) / {length) and we consider the terms of 0{m/L) 
and 0{Gm? /L'^). There will be bilinear contributions from pairs of the (2)— and 
(Q)-components and pairs of the (1)— and (f)-components of the metric. We find 
the former pairs give the spin-induced force and the latter pairs give the radiation 
reaction force. A straightforward computation results in 

— P°(r, r) = 0{G'^rrt'/L^) + (r-dependent terms), (5-39) 

^P\T,r) = -'^R^^,s{z{T))r,{T)i^{T)S^'{T) 

+^e^^,{T)f\{T) (2i"(T)i^(T) +5°^(z(r)) 

+mi°(r)^/,,(T) 

+0(G^m^/L^) + (r-dependent terms). (5-40) 

Taking the T-derivative of Eqs. (5-23) and (5-24), we obtain the equations of motion, 

^z-{T) = -^ep,s{T) {2i^{T)g^^z{T))i\T) - z^{T)i'y{T)g-'{z{T))) 

-\R^p,s{z{T))i^{T)S'^'{T) + 0{G^m^/L^) . (5-41) 

Introducing the proper time r of the orbit, 

^ = 1 + ^0^^(T)i-(T)i'3(T), (5-42) 

we finally arrive at 
D dz"" , , Gra^ , J dz"^ dzf^ dz^ dz^ dz^ ^^ , ,dz^ dz^ dz^ ^. , \ 

--R-^^s{r)^S^\T) + 0{G^myL^) , (5-43) 
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where Q{t) = Q{z{t)). One finds that the result is exactly equal to Eq. (4-54) except 
for the spin-curvature coupling term. 

In the case of a monopole particle discussed in the previous section, the (f) 
matching gave two conditions (4-49) and (4-50). The latter condition was crucial 
to obtain the 0{Gm/L'^) correction terms in the evolution equations of /". In the 
present analysis, we do not have the counterpart of this condition. This indicates 
that the gauge condition relating with the notational mode must be specified to 
determine Dfai/dT. 

§6. Discussion 

Let us first discuss the physical meaning of the equations of motion obtained 
in the preceding two sections. For simplicity, we consider the case of a monopole 
particle. We divide the perturbed metric in the external scheme into the two: 

hfiu{x) = h(^„^^^(,'^^^{x) + h^iaii:)ijLU\X)^ (6T) 

where h(taii)nu{.^) is the part due to the v^^'^^ in the Green function while h(^^ano)fiu 
is due to the Ufj_^ai3 term (see Eq. (3-34)). The singular behavior of the perturbed 
metric in the coincidence limit x ^ z is totally due to hi^jnono)fj.uix)- Thus, we 
introduce the regularized perturbed metric as 

9ireg)fiu{x) ■= 9^lu{x) + h(^taU)fj.u{x), (6-2) 

which has no singular behavior any more. Then we find the equations of motion 
(4-51) and the evolution equations of the triad basis (4-52) coincide with the geodesic 
equation and the geodetic parallel transport equation, respectively, on the regularized 
spacetime with the metric g[reg)fj,u- To see this let us consider the parallel transport 
of a vector A" along a geodesic x" = z"(f) in this spacetime. It is given by 



df - df +^^(reg) P,A ^^ 



^" ■= 3H^" + ^nreg)%'yA''^ = 0, (6-3) 



to the linear order in /i(iaj«)/xi^ where 

^^{reg)°'P'Y '■= ^ [^(tail)"' i3;-f + h(tail)°'r,P ~ ^(tail)fi-t''^j ■ (6-4) 

Then one recovers Eqs. (4-51) and (4-52) by identifying f with T and replacing A"' 
with dz^/dT or /"j. In the case of a spinning particle, there exists an additional force 
in the equations of motion (5-41) due to the coupling of the spin and the background 
curvature. 

The result for the monopole particle seems analogous to that in the electromag- 
netic caseu\ except that the instantaneous reaction force which is proportional to 
higher derivatives of the particle velocity is absent in the present case. This is be- 
cause the particle traces a geodesic in the lowest order approximation. If an external 
force field exists, the assumption of the geodetic motion in the lowest order breaks 
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down and furthermore the contribution of the external force field to the energy mo- 
mentum tensor must be taken into account. Since this fact makes the problem too 
complicated, it is beyond the scope of the present discussion. 

Now let us consider how to construct g{reg)^ivi ™ the case of a monopole particle. 
Unfortunately, we do not have any satisfactory formalism that can be applied to such 
a calculation, even for a specific background spacetime such as the Kerr geometry, 
mainly due to the difficulty in evaluating the bi-tensor v^^a^ix, z). Here we just give 
a few primitive discussions on this matter. 

Basically, there seems to be two approaches for calculating g{reg)p.v (or equiva- 
lently h(taii)nv)- The first one is to calculate h(iaii)nv directly. The second one is to 
calculate h^,y = hi^rnono)^,u + h{tail)^lv and subtract h(^„^ono)iiiy from it. In the following, 
we discuss only the first approach. As for the second approfich, we have nothing to 
mention here, but this direction of research may be fruitful t3^ 

By definition, h(^jnono)iiiy evaluated on the particle trajectory is independent of 
the past history of the particle .til Therefore if we consider the metric defined by 

hlf/\x) = Gm (^6/5,- - ^g,,{x)gf^x)^ j"'^^ dr'G^^t^Cx, z(r'))i"(r')i^(T'), 

(6-5) 
for any finite At {> 0), it will not contain /i(m,ono)^i/ when it is evaluated on the 

particle trajectory. The difference between h\ii/ and h(^faii)nu comes from the integral 
over a small interval, 

r^GmT dr\^^p{x,z{T'))i^iT')z^ir'). (6-6) 



Since Vpaap{x-,z) is regular in the coincidence limit x ^^ z, this integral will be 

negligible for a sufficiently small At. Thus lim^T-^o h\ii,^ will give h(^iaii)iiu- 

In the case of the electromagnetic (vector) Green function, a calculation along 
the above strategy was performed by DeWitt and DeWittlHI^ by assuming the back- 
ground gravitational field is weak so that its metric is given by the small perturbation 
on the Minkowski metric, 

gfiu = Vfj-iy + ^IS • (6-7) 

DeWitt and DeWitt calculated the relevant part of the Green function perturbatively 
to the first order in h^il by using the Minkowski Green function. 

Here we should mention one important fact. We have obtained the equations of 
motion with the correction term of 0{Gm/L'^). Although we use the terminology 
'radiation reaction' to describe it, it is not appropriate in a narrow sense because 
the correction term may well contain something more than just the usual effect of 
radiation reaction. In fact, in the electromagnetic case, the existence of the effect 
which is termed as 'the induced polarization force on the background spacetime' is 
reported by several authors t3'. Furthermore, a calculation analogous to that done 



*' There is a possibility that the future light cone emanating from z crosses the particle trajectory 
again. Since inclusion of this possibihty makes the problem too complicated, we do not consider it 
here. 
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by DeWitt and DeWitttll) for the electromagnetic case was done by Carmelill^ for 
the gravitational case and it was shown that the tail part correctly reproduces the 
lowest order post-Newtonian corrections to the equations of motion. However, no 
such calculation has been done for the background with strong gravity, such as a 
black hole spacetime. It seems difficult to develop DeWitt and DeWitt's approach 
to higher orders in h^il. It is a challenging issue to formulate a systematic method 
to evaluate the tail part of the metric when the background gravity is strong and 
clarify its physical content. 

Turning back to the effect of the gravitational radiation reaction, we should make 
one additional comment. There has been some proposals to obtain the radiation 
reaction force in a quite different manner. Among others is the use of the radiative 
Green function (a half of the retarded minus advanced Green functions) in the case 
of a Kerr background proposed by Gal'tsovlllP. As easily seen from the results in 
section |3|, the use of the radiative Green function instead of the retarded one results 
in the replacement of ip(taii)fiu{^) by '^^JwC^)' which is defined by 

/ + 00 
dT'vp,^,p,{x,z{T'))i»{T')i''{T'). (6-8) 

Gal'tsov proved that the back reaction force computed using the radiative Green 
function gives the loss rates of the energy and the z-component of the angular mo- 
mentum of the particle in quasi-periodic orbits which correctly balance with the 
emission rates of the corresponding quantities by gravitational radiation. However, 
we do not think that this fact indicates the correctness of the prescription, even if we 
restrict it to the case of a Kerr background, because those constants of motion are 
special ones which reflect the existence of the corresponding Killing vector fields. For 
such quantities, there may be some symmetry in the structure of the Green function 
which makes the use of the radiative Green function valid. However, it is doubtful 
that the radiative Green function correctly describes the radiation reaction effect on 
the Carter constant. 

Finally we make a couple of comments on the implications of our results. It 
is important to note that the particle does not have to be a black hole since the 
detailed internal boundary condition was not used to determine the metric in the 
internal scheme. The resulting equations of motion should be equally applicable to 
any compact body such as a neutron star. The essential assumption here is that the 
only length scale associated with the particle is Gm. In this sense, we have shown 
the strong equivalence principle to the first order in Gm. 

We also note that our results strongly support, if not rigorously justify, the so- 
called black hole perturbation approach. In the black hole perturbation approach, 
one calculates the gravitational radiation from a particle orbiting a black hole with 
the assumption that the particle is a point-like object with the energy momentum 
tensor described by the delta function. Although this approach has been fruitful, 
there has been always skepticism about the validity of the delta functional source. 
What we have shown in this chapter is the consistency of using the delta function 
in the source energy momentum tensor within the order of matched asymptotic 
expansion we have examined. 
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